Mon. Not. R. Astron. Soc. 000, ??-?? (2009) Printed 29 June 2009 (MN WI^ style file v2.2) 



Too large and overlooked? Extended free-free emission 
towards massive star formation regions 

S. N. M. G. Burton^, E. Keto\ S. Kurtz^ & A. J. Walsh^ 

^Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 
^School of Physics, University of New South Wales, 2052, Sydney, Australia 
^Australia Telescope National Facility, CSIRO, Epping 1710, Sydney, Australia 

^CRyA, Universidad Nacional Autonoma de Mexico, Apartado Postal 3-72, 58089, Morelia, Michoacdn, Mexico 
^ Centre for Astronomy, James Cook University, Townsville, QLD 4811 Australia 



in prep. 



ABSTRACT 

We present Australia Telescope Compact Array observations towards 6 massive 
star formation regions which, from their strong 24 GHz continuum emission but no 
compact 8 GHz continuum emission, appeared good candidates for hyper-compact Hii 
regions. However, the properties of the ionised gas derived from the 19 to 93 GHz 
continuum emission and HTOa + H57a radio recombination line data show the major- 
ity of these sources are, in fact, regions of spatially-extended, optically-thin free-free 
emission. These extended sources were missed in the previous 8 GHz observations due 
to a combination of spatial-filtering, poor surface brightness sensitivity and primary 
beam attenuation. 

We consider the implications that a significant number of these extended Hii 
regions may have been missed by previous surveys of massive star formation regions. 
If the original sample of 21 sources is representative of the population as a whole, the 
fact that 6 contain previously undetected extended free-free emission suggests a large 
number of regions have been mis-classified. Rather than being very young objects 
prior to UCHii region formation, they are, in fact, associated with extended Hii 
regions and thus significantly older. In addition, inadvertently ignoring a potentially 
substantial flux contribution (up to ^^O.SJy) from free-free emission has implications 
for dust masses derived from sub-mm flux densities. The large spatial scales probed 
by single-dish telescopes, which do not suffer from spatial filtering, are particularly 
susceptible and dust masses may be overestimated by up to a factor of ~2. 

Key words: stars: early type, stars: formation, stars: evolution, ISM: evolution, 
(ISM:) HII regions, radio continuum: ISM, radio fines: ISM 



1 INTRODUCTION 

A definition often used to describe a star as 'massive' re- 
quires it to have sufficient UV radiation to ionise the sur- 
rounding environment. The ionised material can then be ob- 
served as free-free continuum emission at radio wavelengths. 
When a massive star is forming, the emergence and subse- 
quent expansion of detectable free-free emission from the 
natal cocoon, provides an indicator of the stars evolution- 
ary state. As the ionising radiation overcomes quenching due 
to the large mass of infalling material, the continuum emis- 
sion emerges. At first the ionised gas is extremely dense, 
but once the ionised gas pressure is sufficient to overcome 



the gravitational potential, the Hii region expands into ul- 
tra compact (UC), compact and finally diffuse Hii regions. 
The spectral energy distribution (SED) at cm wavelengths 
falls into two distinct regimes with increasing frequency. As 
the free-free emission turns from optically thick to thin, the 
spectral index (S^ cx i/" , where Si/ is the flux density at 
frequency i^) changes from a ~ 2 to -0.1. The frequency at 
which the turno ver occurs sc ales almost linearly with the 
electron density (iKurt j [20051 ) i.e., denser and presumably 
younger objects will exhibit the turn over at higher frequen- 
cies. At shorter wavelengths the SED quickly becomes dom- 
inated by thermal emission from dust in the surrounding 
natal cocoon, which has a spectral index, a ~ 3 — 4. 
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To identify sources at the ear liest stages of massi ve star 
formation we used the work of IWalsh et al.l l| 19981 . here- 
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after WBHR98) to select objects in the relatively unex- 
plored southern sky with: (i) SED's peaking in the far- 
infrared, (ii) sufficient luminosity to expect stars massive 
enough to ionise the surrounding material, (iii) no detectable 
cm free-f ree emission and (iv) class II methanol maser 
emission (jLongmore et al ] l2007al . hereafter L07). Class II 
methanol masers are generally found offset from cm contin- 
uum emission but always coincident with sub-mm cores (e.g. 
I Walsh et al.ll2003l ) which have sufficient luminosity to con- 
tain massive stars (WBHR98). These highly luminous cores 
with no detected cm emission are interpreted as massive star 
forming cores prior to the onset of UCHii formation. Subse- 
quent multi-wavelength studies have since derived the large- 
scale properties (mass, chemical composition, temperature 
etc.) of these star forr ning cores using the SEST, Mopra 
and ATCA telescopes /Hill et al j | 2005riPurceU et alllioo^ . 
I2OO9I : iLongmore et al.| [2007b. 200i) and performed a census 
of the (proto-) stellar popul ations at infrared wavelengths 



Table 1. Pointing centre from the methanol maser position, 
Vlsr fc IRAS name (fro m WBHR98) and the adopted distance 
(from iPurcell et alll2006l) of the observed regions. 



or tne (proto-) stellar popui a 
jLongmore et al.ll2006l . boOTj ) . 



In L07 we selected 21 regions for further study, all 
of which contained methanol maser emission detected by 
WBHR98. Six of these regions were reported as containing 
continuum emission by the simultaneously-observed 8 GHz 
observations of WHBR98, while the remaining 15 were re- 
ported as devoid of 8 GHz continuum emission by their 
8 GHz observations. Using this selection criteria, we hoped 
to identify 15 young regions prior to UCHII formation and 
6 older regions, post-UCHII formation. However, contrary 
to expectation, twelve rather than six of these regions were 
detected in the 24 GHz continuum by L07: the six detected 
in the 8 GHz continuum by WBHR98 and six reported as 
devoid of 8 GHz continuum by WBHR98. 

In this work we are seeking to find the nature of the 
sources detected at 24 GHz by L07 but not at 8 GHz by 
WBHR98. If these cores are indeed at an earlier evolution- 
ary stage as originally thought, this emission could be free- 
free emission in the optically-thick regime between 8 and 
24 GHz. In this case their 24 GHz flux should be ~10 times 
brighter than at 8 GHz and may explain the 8 GHz non- 
detections. However, both L07 and WBHR98 observations 
were made using an interferometer with significantly dif- 
ferent array configurations and so probed the source struc- 
ture at different spatial frequencies and with different sur- 
face brightness sensitivities. The other possibility discussed 
in L07 is that these sources could instead be much older, 
spatially extended, optically thin Hii regions which were 
resolved out/not-detected d ue to the extended array con- 
figuration used by WBHR98. feurtz et all ()l999[ ) found just 
such extended continuum emission surrounding 12 out of 15 
UCHii regions they observed, with important consequences 
for interpreting the nature of these regions. If such a result 
was found for many of the L07 sources, and these were repre- 
sentative of the WBHR98 sample as a whole, this would have 
ramifications regarding their currently assumed age and evo- 
lutionary state. Not only would this imply the sources are 
significantly older than originally thought, but it may also 
affect attempts to derive dust masses at shorter wavelengths 
ignoring the contribution from the extended free-free emis- 
sion. In this paper we aim to uncover the nature of these 
sources through continuum and radio recombination line ob- 
servations from 19 to 93GHz. 
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2 OBSERVATIONS AND DATA REDUCTION 

The sources were selected from L07 as those with 24 GHz 
continuum detections but no corresponding 8 GHz detec- 
tions reported in WBHR98. Table [T] lists the pointing cen- 
tres from the m ethanol maser posit ion, Vlsr and kinematic 
distances (from [Purcell et al.ll2006l ) for each of the regions. 



2.1 ATCA continuum and radio recombination 
line (RRL) data 

The Australia Telescope Compact Array (ATCA) observa- 
tions were undertaken between 2006 and 2008 as part of 
ATCA projects C1287 and C1751. The observing setups 
(array and correlator settings, approximate sensitivities and 
beam sizes, continuum frequencies and RRL transitions) are 
listed in Table [2] In order to mitigate problems constructing 
SEDs of extended sources with different beam sizes, three 
array configurations were used (H75, H168 and H214, with 
maximum baselines of 75m, 168m and 214m, respectively) 
in an attempt to match the beams at the different observing 
frequencies. These configurations contain both north-south 
plus east-west baselines to allow for snapshot imaging. Each 
source was observed in multiple 15 minute cuts, with a close 
phase calibrator observed for 3 minutes before and after 
each cut. PKS 1253-055 and PKS 1921-293 were used as 
the bandpass calibrators. The ATCA standard PKS 1934- 
634 was used as the primary calibrator for observations be- 
low 30 GHz and Uranus was used at higher frequencies. The 
data were reduced in the standard way using the MIRIAD 
(jSault et al.lll995l ) package. Bad visibilities were fiagged and 
the gains/bandpass solutions from the appropriate calibra- 
tor were applied to the visibilities. A zeroeth order polyno- 
mial was fit to the line-free channels and subtracted from 
the visibilities to separate the line and continuum emission. 
These separate visibilities were then Fourier transformed to 
form image cubes. Due to variation in system temperatures 
between the antennae, individual visibilities were weighted 
by their system temperature. The images were cleaned and 
convolved with a Gaussian synthesised beam to produce a 
final restored image. 

Continuum fiuxes were extracted using imfit on the re- 
stored images. Extended sources were fit with a 2D Gaus- 
sian profile to calculate the spatial extent, integrated and 
peak fiuxes densities. Unresolved emission was fit as a point 
source. The positions, peak fiux and deconvolved sizes of the 
continuum emission from these fits at each frequency and 
array configuration are shown in Table [S] Where detected. 
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RRL spectra were extracted for each source at the peak of 
the continuum emission, baseline subtracted and fit with a 
Gaussian profile using the gauss method in CLAS£|. Ta- 
ble |4] lists the velocity-integrated fiux, Vlsr, FWIfM, peak 
flux and their corresponding uncertainties for each of the 
detected transitions and array configurations. 



3 RESULTS 

§ 13.11 to 13.61 report the 19 to 93 GHz continuum and RRL 
detections towards each of the individual regions. These are 
compared to both the L07 24 GHz observations and original 
8 GHz WBHR98 data to try and uncover the nature of the 
continuum sources detected by L07 but not WBHR98. We 
focus on analysis of the continuum rather than the RRL 
data, which will appear in a subsequent paper. In this work 
we simply refer to the RRL detections as further evidence 
that the continuum arises from free-free emission. 

3.1 G316.81-0.06 

3.1.1 Previous observations 

As shown in Figure [T] (top-left), L07 detected three sources 
in the region; two 24 GHz continuum sources (core 2 at 
the methanol maser position, and core 3 offset ~30" south- 
west from the maser) and an ammonia filament extending 
east from core 2 and the methanol maser position. Despite 
the strong (>0.5Jy) continuum towards core 2 at 24 GHz, 
WBHR98 only reported a single 8 GHz detection towards 
this region at the position of core 3. To investigate this 
apparent discrepancy between the L07 and WBHR98 ob- 
servations, we first re-examined the WBHR98 8 GHz data. 
With poor uv-coverage, phase stability and strong sidelobe 
emission, the image fidelity was not sufficient to rule out an 
8 GHz detection at the position of core 2. However, the si- 
multaneously observed WBHR98 6.7 GHz continuum images 
have much better image fidelity and show a non-detection 
to 5 mjy at the position of core 2 and a 40 mjy detection at 
core 3 with a similar morphology to the 24 GHz detection. 
The apparent discrepancy in the measured flux density of 
core 2 between the lower-frequency, very extended array ob- 
servations of WBHR98 and the higher-frequency, compact 
observations of L07 therefore appears to be related to the 
nature of the source emission. Based on the new observa- 
tions, we investigate the nature of core 2 below. 

3.1.2 G316.81-0.06: Core 2 

As shown in Figure [TJ we confirm the strong 24 GHz contin- 
uum L07 detection of G316.81-0.06 core 2. At higher reso- 
lution the cm-continuum emission is resolved into two com- 
ponents separated by ~10" north-south. At 93 GHz we find 
unresolved continuum emission at the centre of the two cm- 
continuum components. 

It is clear from the 19, 24 and 35 GHz continuum im- 
ages (Figure [TJ and flux densities (Table [3| , that the bulk 
of the emission in core 2 is spatially extended. Although 

1 http://www.iram.fr/GILDAS 



the more extended H168 array configuration 35 GHz obser- 
vations resolve the north-south elongation into at least two 
components, these only make up a small fraction of the fiux 
compared to that measured at the same frequency in the 
more compact H75 array configuration. 

Strong emission from both H70Qf and H57a are detected 
towards core 2 in all the array configurations (Figure [2]). 
Similar to the cm-continuum emission, the peak flux density 
of the H57q: RRL in the more extended array configuration 
only recovers a small fraction of the flux density in the more 
compact conflguration. 

To correct for this spatial filtering when creating an 
SED it is therefore necessary to sample the source emission 
at the same spatial frequencies. Figure [13] shows the SED 
for core 2 sampled between spatial frequencies of 5— lOkA, 
the mutually overlapping range between all the 18.9 to 
35 GHz observations. This shows that on spatial scales of 
25-50", the spectral index is almost flat - consistent with 
optically-thin free-free emission. It is not possible to unam- 
biguously compare this spectral index with the flux from the 
92 & 93 GHz observations due to the smaller primary beam, 
much higher spatial frequencies and the potential contri- 
bution from thermal dust emission. However, the similarity 
between the 92 GHz H75 array conflguration observations 
and the lower frequency emission suggests at least some of 
the 92 GHz emission may be due to free-free emission. 

The 93 GHz observations using the more extended H214 
array configuration show a very different structure, with only 
a single, weak, unresolved detection directly in the middle of 
the two compact lower frequency knots. This source is also 
coincident with an N2H"'" core (Longmore et al. in prep) trac- 
ing dense molecular gas, so we postulate this 93GHz emis- 
sion is likely to be from thermal dust. Although the emission 
from this source is probably distinct from that of core 2 at 
lower frequencies, they are spatially associated within the 
beam of the lower resolution data, and as such are both 
labelled as 'core 2'. 

Finally, there is also a strong 93 GHz continuum source 
at the position of the ammonia peak, core 1. The spectral 
index inferred from the upper limits of the lower frequency 
data at this position means this is also likely to be thermal 
dust emission. As this source is substantially outside of the 
primary beam, the flux measurement reported in Table [3] is 
highly uncertain. 

3.2 G331.28-0.19 

Figure [3] shows images of the continuum emission detected 
towards the region. The morphology of the 19 and 35 GHz 
continuum emission is very similar to that observed at 
24 GHz by L07: a single source is seen peaking close to the 
methanol maser position with an extended tail towards the 
east. This continuum emission is offset to the south east 
from an L07 ammonia core which is extended north-east to 
south-west. 

WBHR98 reported no 8 GHz continuum source in the 
region. Re-examination of the WBHR98 8 GHz data shows 
a potential source of 4mJy at this position with a simi- 
lar morphology. However, due to poor weather, only 2x2 
minute cuts were observed in the WBHR98 observations so 
the image fldelity is poor and the emiss ion brightness l ay be - 
low their detection criteria. However, [Phillips et al.l (|l998l ) 
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Table 2. Overview of the observing setup. For the H70«, H57« and H41« observations, the observing frequencies were centred on the 
respective RRL transitions and continuum emission was extracted from Une-free channels. The 93.5 GHz observations were solely high 
spatial resolution continuum observations and no RRL was observed in this setup. Columns three and four give the bandwidth (BW) 
and number of channels for each of the correlator settings used. AV shows the corresponding velocity resolution per channel for each IF 
in kms~^. Column six lists the array configuration for each observation. The H75, H168 and H214 array configurations have baselines 
ranging from 31 - 89 m, 61 - 192 m and 82 - 247 m, respectively. Columns 7 and 8 show the theoretical spectral line RMS per channel and 
theoretical continuum brightness sensitivity (baaed on a one hour observation in average weather conditions using the ATCA sensitivity 
calculator - http://www.atnf.csiro.au). The final column gives a characteristic synthesised beam size (for a source at 5 = —30°). 
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report weak (~3.5mJy) continuum emission at 8 GHz, ex- 
tended east- west, with a similar morphology to that in Fig- 
ure [S] 

The H70q: and H57Qf RRL emission detected towards 
the continuum (Figure |4]) and spectral index of —0.07 de- 
rived from the SED (Figure [13} confirms the continuum be- 
tween 19 and 35 GHz as free-free emission. 

A compact continuum source is also detected at 92 and 
93 GHz. This is offset from the peak of the lower frequency 
emission and is instead associated with the methanol maser 
emission, the peak of the L07 NH3(1,1) ammonia emission 
and the unresolved NH3(4,4) and (5,5) emission. Due to a 
combination of the lower sensitivity to extended emission of 
these observations and the lower optical depth of free-free 
emission at higher frequencies, the 92 & 93 GHz emission 
could be from the densest and brightest free-free component 
inside the larger region traced by the lower frequency emis- 
sion. However, given the similar synthesised beam sizes of 
the 35 GHz H168 and 92 GHz H75 observations, it is diffi- 
cult to understand the offset in the peak of the emission. We 
favour an alternative explanation in which the 92 & 93 GHz 
emission is dominated by thermal dust emission rather than 
free- free emission. 

In summary, G331.28-0.19 core 1 appears to be an opti- 
cally thin HH region. The 92 GHz emission and the ammonia 
emission are likely tracing a hot molecular core located at 
the edge of the HH region, similar to the case of G29.96- 
0.02 (|Cesaroni et al.lll998l ). The 8 GHz non-detection by 
WBHR98 can be attributed to their poorer sensitivity for 
this source compared to the rest of their sample. 

3.3 G12.68-0.18 

Figure [5] shows images of the continuum emission detected 
towards the region. The morphology of the 19 and 35 GHz 
continuum emission is similar to that observed at 24 GHz 
by LOT: a single extended source is seen 10" north-east of 
the methanol maser emission and offset from the ammonia 
cores in the region. Figure |S] shows the H57q: and H70q: 
RRL spectra detected towards the peak of the continuum 
emission, confirming this arises from ionised gas. Figure [T3l 
shows the SED for this source. Although the quality of the 
fit is not particularly good, the spectral index of —0.4 ±0.11 
is clearly inconsistent with optically thick free-free or dust 



emission. The extremely negative spectral index is biased by 
the low fiux density of the 35 GHz H168 data point, which 
is probably due to missing fiux. Removing this point before 
fitting results in a spectral index of —0.18 ± 0.09. 

Re-inspection of the WBHR98 data shows this source is 
~10' from the WBHR98 pointing centre, significantly out- 
side their half-power beamwidth of 5.7'. Any emission would 
therefore have suffered from primary beam attenuation in 
their observations. 

Figure [S] also shows a second, unresolved source de- 
tected at 93 GHz, 5" south of the maser position, coincident 
with the unresolved NH3(1,1) (5,5) core reported in L07 
(core 4). The non-detection at frequencies <93 GHz makes 
it likely that the emission mechanism for this source is from 
dust rather than thermal free-free emission. 

In summary, G12.68-0.18 core 3 is an extended region of 
optically-thin ionised gas that was not detected by WBHR98 
as it lies far outside their primary beam. 

3.4 G19.47+0.17 

Figure [7] shows images of the continuum emission detected 
towards the region. The morphology of the 19 and 35 GHz 
continuum emission is similar to that observed at 24 GHz 
by L07: an unresolved continuum source, labelled core 2, 
lies ~15" north-west of the methanol maser position. H57Qf 
and H70a RRL emission is detected towards the peak of 
the continuum emission (Figure |8]), confirming the emission 
arises from ionised gas. The spectral index (Figure I13|l is 
consistent with optically-thin free-free emission. 

No 8 GHz continuum emission was reported at this lo- 
cation by WBHR98. Re-examination of the WBHR98 data 
shows this region lies outside the primary beam so the con- 
tinuum sensitivity would have suffered significantly from pri- 
mary beam attenuation. 

A second, weaker source is seen only at 92 and 93 GHz 
within a few arc-seconds of the methanol maser position and 
unresolved NH3(1,1) — > (5,5) emission, designated core 1 in 
L07. At 35 GHz the 3a upper limit is 7 mjy. This implies a 
spectral index between 35 and 92 GHz of a > 2, which rules 
out either optically-thick or thin bremsstrahlung emission. 
We conclude that the core 1 emission is most likely from 
warm dust. 
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Table 3. Continuum detections towards each of the regions. The frequency and array configurations are given in columns two and 
three. The Usted positions, fluxes and angular sizes are taken from the fits to the images, as outlined in §|2] For extended sources, the 
angular sizes have been deconvolved with the synthesised beam. The errors shown in column seven are from the fit to the data and 
do not take into account the absolute error in the flux scale which is assumed to be 20% and 30% for the observations < 90 GHz and 
> 90 GHz, respectively. " the integrated flux for extended sources is shown in column seven in parentheses. * for integrated fluxes (in 
parentheses) the units are mjy rather than mjy/beam. G316.81-0.06 Core 1 is considerably outside the primary beam of the 92 GHz 
H75 observations making this flux measurement highly uncertain. 
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Table 4. Detected radio recombination line emission towards each region. The second and third columns give the transition and array 
configuration. Columns four through seven list the velocity-integrated line flux, central velocity, full-width at half maximum and peak 
intensity with associated errors from fits to the spectra (see §[2]l, extracted from the position of peak continuum flux in Table |3] Column 
8 gives the line to continuum ratio for each transition. " the reported flux is the velocity-integrated flux extracted from the peak spectra. 
The spectra are shown for each region separately in Figures |21 l4l [6l [8l 1101 and 1121 and overlayed with the fit to profile. 



Source 


Trans. 


Array 


Flux" 


Vlsr 


FWHM 


Peak 


Linc/cont. 








(mjy/beam.km/s) 


(km/s) 


(km/s) 


(mjy/beam) 


ratio 




no 1 (X 


117^ 




-^o.uinu.uo 




ouu 


u.o 




H57a 


H168 


1660±54 


-46.5±0.34 


20.8±0.84 


72 


0.5 




H70o 


H168 


2483±26 


-43.8±0.11 


20.8±0.26 


110 


0.3 


G316.81-0.06 Core 3 


H57o 


H75 


3695±62 


-43.2±0.20 


24.1±0.49 


140 


0.5 




H57o 


H168 


1630±48 


-45.5±0.28 


19.4±0.68 


79 


0.5 




H70a 


H168 


2737±32 


-42.3±0.14 


24.8±0.35 


100 


0.2 


G331.28-0.19 Core 1 


H57o 


H75 


1363±19 


-76.6±0.12 


17.9±0.28 


71 


0.5 




H70o 


H168 


429±20 


-78.6±0.48 


21.8±1.2 


18 


0.3 


G12.68-0.18 Core 3 


H57ci 


H75 


450±23 


57.5±0.85 


33.8±2.0 


13 


0.4 




H57o 


H168 


273±22 


55.8±0.98 


22.7±2.0 


11 


0.5 




H70o 


H168 


780±64 


54.7±1.3 


34.3±3.6 


21 


0.3 


G19.47-I-0.17 Core 2 


H57a 


H75 


429±53 


16.4±1.2 


20.9±3.2 


19 


0.6 




H70o 


H168 


210±37 


12.7±1.4 


15.4±3.0 


13 


0.2 


G24.79+0.08 Core 1 


H57a 


H75 


1235±80 


113.7±0.80 


26.2±2.1 


44 


0.4 




H57o 


H168 


735±71 


112.8±1.3 


27.2±3.2 


25 


0.4 


G24.79-I-0.08 Core 3 


H57o 


H75 


800±62 


109.9±0.58 


16.4±1.6 


46 


0.6 




H70ci 


H168 


569±30 


112.4±0.45 


17.3±1.1 


31 


0.2 


G24.85-I-0.09 Core 1 


H57o 


H75 


560±60 


111.9±1.2 


24.1±3.2 


22 


0.4 



3.5 G24.79+0.08 

Figure |5] shows images of the continuum emission detected 
towards the region. The morphology at 19 and 35 GHz is 
very similar to that observed at 24 GHz by L07; an unre- 
solved source at the methanol maser position (core 1) and 
a slightly extended source at the edge of the primary beam 
(core 3). Although, core 3 was not detected at 92 & 93 GHz, 
it lies outside the primary beam at >24 GHz making both 
detection more difflcult and the flux determination more un- 
certain. Figure [TUl shows H57a and H70a RRL emission is 
detected at the peak of the continuum emission towards core 
3 and H57a is detected towards core 1, confirming both these 
sources contain ionised gas. 

WBHR98 reported a la non-detecti on to 3mJy at both 
24 GH z continuum positions. However, iForster fc Caswelil 
reported two 8.6 GHz continuum sources in the prox- 
imity of the methanol maser emission, with flux densities of 
36 and 51 mjy. Re-examining the calibrated WBHR98 uv 
data shows residual phase errors in the data for this region. 
The image shows 8 GHz continuum emission in the region 
but its location is unclear, resulting in the reported non- 
detection. 

Core 1 is well know n in the literature as a hyper- 
compact Hii region (e.g. iBeltran et all |2007| ). Figure 1131 
shows the SED at the frequencies observed in this work 
which are consistent with the source transitioning be- 
tween optically-thick and thin free-free emission. The single- 
component spectral index attributed to this source in Fig- 
ure [13] is therefore not a good representative of the true 
SED. 

In summary, these observations confirm both the L07 
24 GHz continuum detections and show core 3 is an ex- 
tended region of optically-thin ionised gas which was missed 
by WBHR98. Core 1 is a well known optically-thick, hyper- 



compact Hii region which was missed by WBHR98 due a 
combination of the steeply rising fiux density between 8 and 
24 GHz and calibration problems. 

3.6 G24.85+0.09 

Figure [TT] shows images of the continuum emission detected 
towards the region. The morphology at 19 and 35 GHz is 
similar to that observed at 24 GHz by L07: a compact source 
is seen peaking ~5" west of the methanol maser position, 
with a fainter tail extending towards the south-east. No 
emission was detected at 92 or 93 GHz. The detection of 
H57a (Figure [T^ confirms the presence of free-free emission 
towards the source. Figure [T51 shows the SED. The spectral 
index of 0.4 determined from the fit is not very robust due 
to the flux uncertainties and the narrow range of detected 
frequencies. Nevertheless, the spectral index is sufficient to 
rule out dust as the primary contributor to the emission. The 
spectral index of 0.4 is intermediate between the canonical 
values of -\-2 and —0.1, suggesting that the sampled frequen- 
cies are near the turn over region. Re-examination of the 
WBHR98 data shows a 3(j non-detection to 6mJy at 8 GHz. 
Extrapolating the WBHR98 upper limit assuming a spectral 
index of -1-2 results in a 19 GHz flux density well below the 
measured value (see Figure [T3}. From this, and the extended 
source size we assume the emission is optically thin. 

3.7 Spectral Energy Distributions 

Figure [13] shows the SEDs towards each of the cores de- 
tected at 24 GHz by L07 and not at 8 GHz by WBHR98. 
Our aim was to distinguish between the different expected 
emission mechanisms: optically-thin/thick thermal free- free 
emission (a ^ —0.1 or 2, respectively. Given the different 
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Table 5. Nature of the L07A and newly detected continuum 
sources towards each of the regions derived from the spectral en- 
ergy distributions (see i|3.7l l. 



Region 


Core 


Nature of emission 


Freq. 


G316.81-0.06 


2 


optically thin free-free 


24 GHz 


G331. 28-0.19 


1 


optically thin free-free 


24 GHz 


G12.68-0.18 


3 


optically thin free-free 


24 GHz 


G12.68-0.18 


4 


dust 


92 GHz 


G19.47+0.17 


1 


dust 


92 GHz 


G19.47+0.17 


2 


optically thin free-free 


24 GHz 


G24.79+0.08 


1 


optically thick free-free 


24 GHz 


G24.85+0.09 


1 


optically thin free-free 


24 GHz 



range of spatial frequencies probed between the different 
observations, the observations at <90 GIfz were re-imaged 
selecting only visibilities in the mutually overlapping spatial 
frequency range. In practice this meant dropping the longest 
baselines at higher frequency, thus biasing the SED towards 
the large scale structure. As a result it was not possible to 
resolve G316.81-0.06 Core 2 into the two separate compo- 
nents, and only a composite SED is given. After re-imaging 
and re-fitting the emission to determine the flux density (as 
outlined in §[2]), the spectral index was calculated using a 
weighted, linear least-squares fit. Uncertainties in the indi- 
vidual flux density measurements were calculated by adding 
in quadrature the formal error of the flt to the emission and 
the estimated uncertainty in absolute flux (see § [2]). The 
absolute flux uncertainty was dominant in nearly all cases. 
As there were no overlapping spatial frequencies between 
the higher resolution observations at 92/93.5 GHz and the 
remaining data, it was not possible to re-image these data 
in the same way. As a result, the 92/93.5 GHz observations 
miss much of the flux seen towards extended sources at lower 
frequencies, potentially biasing the spectral index. The un- 
certainty to the spectral index is provided in each of the 
figures. 

Table [5] shows the conclusions as to the nature of the 
24 GHz continuum emission in L07 based on the SED fit- 
ting. Despite the sometimes large uncertainty in spectral in- 
dex, the SEDs are sufficient to rule out an optically-thick 
component as the dominant emission mechanism at this 
spatial scale in G316. 81-0.02 core 2, G331.28-0.19 core 1, 
G12.68-0.18 core 3 and G19.47-I-0.17 core 2. This does not 
preclude the existence of compact ionised gas components 
at smaller spatial scales (indeed, density gradients may well 
exist), only that at this position and spatial scale diffuse, 
optically thin ionised gas is dominant. 

So how do these SEDs compare to the WBHR98 data? 
Where applicable, the WBHR98 8 GHz upper limits are 
shown as arrows in Figure 1131 Given the sometimes sev- 
eral orders of magnitude discrepancy between the expected 
8 GHz fiux extrapolated from the SEDs, the non-detections 
seem initially inconsistent. One potential explanation may 
be that the emission is optically-thick between 8 and 19 GHz, 
and as a result drops below the WBHR98 detection limit. 
To test this we have extrapolated the expected SED be- 
tween 8 and 19 GHz assuming the emission is optically thick 
and just falls below the WBHR98 detection limit (dashed 
line in Figure I13|) . While this may be a viable explana- 
tion for G24.79-f0.08 core 1 and possibly G24.85+0.09, it 



clearly fails for the remaining sources. We conclude the only 
other viable explanation is that WBHR98 observations must 
have missed the spatially extended optically-thin emission 
through either spatial filtering, reduced surface brightness 
sensitivity or a combination of the two. 



4 DISCUSSION: THE "MISSING" 
CM-CONTINUUM SOURCES 

Observations of cm-continuum emission towards regions of 
massive star formation provide a powerful indicator of the 
region's age. Simplistically, detection of relatively weak (1- 

5 mjy) emission from shock-ionised gas signals the existence 
of a powering source (that may be driving a thermal jet 
for example), while stronger emission from photo-ionised 
gas signals a young massive stellar object has evolved suffi- 
ciently to have begun ionising its surrounding environment. 
Conversely, lack of bright cm-continuum emission is often 
used to infer the region is at an evolutionary stage prior 
to the development of an UCHii region. As such, previous 
large, high-resolution, cm-continuum surveys have proven 
invaluable in helping to assess the relative ages of many 
massive star formation regions. However, from our original 
sample of 21 massive star formation regions, selected from 
such a survey to contain no cm-continuum emission, we have 
confirmed that 6 of these are in fact associated with cm- 
continuum emission. In 14.11 we summarise the reasons for 
these discrepancies for the sources we have studied in this 
work and in § 14.21 discuss the implications this has on the 
implied evolutionary stage of the objects. We then assess 
how these results may affect the interpretation of the whole 
WBHR98 sample. 



4.1 Reasons for the "missing' cm-continuum 
sources 

We find three reasons why the cm-continuum sources in this 
work were not detected at 8 GHz by WBHR98. These can 
be separated by the nature of the source emission and the 
sensitivity limit of the survey as follows: 

(i) the free-free emission is dominated by compact com- 
ponents, is optically-thick at 8 GHz and lies below the 
WBHR98 detection limit (i.e. G24.79-I-0.08 core 1) 

(ii) the emission is significantly above the WBHR98 
detection limit but much more spatially extended than 
the WBHR98 beam (i.e. G316. 81-0.06, G331.28-0.19 and 
G24.85-f0.09) 

(iii) the emission lies well outside the primary beam of 
WBHR98 (i.e. G12.68-0.18 and G19.47-h0.17) 

In each of the three cases, the reason for the non- 
detection at 8 GHz is straightforward: (i) is due to the sen- 
sitivity hmits of WBHR98; (ii) is due to the WBHR98 ob- 
servations resolving-out extended emission, and, (iii) is due 
to the significantly reduced sensitivity of the observations 
outside of the primary beam. 
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4.2 Implications 

4-2.1 Methanol masers and an evolutionary sequence 

Based on the WBHR98 data, 15 of the 21 regions in the ini- 
tial L07 survey were selected as being methanol maser sites 
devoid of 8 GHz continuum emission and were thus assumed 
to be in an evolutionary stage prior to UCHii region forma- 
tion. The results from this work show that in fact only 9 of 
the original 21 regions contain no free-free emission and the 
remaining 12 are either at the HCHii region stage [case (i)] 
or in fact substantially older [cases (ii or iii)]. 

WBHR98 found ~80% of their methanol maser sites 
were not associated with radio continuum emission and sug- 
gested methanol masers may form before an UCHii region 
develops. While this work does not contradict the conclu- 
sions of WBHR98 - methanol masers do trace very early 
stages of massive star formation - our results suggest a 
substantially larger fraction of methanol masers than de- 
termined by WBHR98 are associated with more evolved 
stages, UCHii regions and later. However, it is not possi- 
ble to rule out methanol masers as signposts of the pre-HII 
phase. Given the clustered mode of massive star formation, 
these youngest sources may simply be coincident with more 
evolved massive young stellar objects. 

Our findings support the idea that methanol masers 
commonly occur in both a very early phase before a de- 
tectable UCHii region develops (this may include the hy- 
percompact Hii region phase), and after the formation of 
a UCHii region where the ionised gas is expanding and 
becomes extended enough that it may be missed by ex- 
tended interferometer configurations, such as that used by 
WBHR98, due to spatial filtering/poor surface brightness 
sensitivity. 

In order to fully understand the relative occurrence of 
methanol maser sites in these three stages, we suggest a 
radio continuum survey that is sensitive to more extended 
radio continuum emission be undertaken. Such a survey 
should observe a large number of methanol maser sites with- 
out associated radio continuum emission, as identified by 
WBHR98. 

4-2.2 Dust masses derived from sub-mm flux densities 

Inadvertently ignoring a potentially substantial flux contri- 
bution (up to ~0.5Jy) from free-free emission has implica- 
tions for deriving dust masses from sub-mm flux densities. 
Given the sharply rising flux density of dust with frequency 
this overestimate will be most pronounced at lower frequen- 
cies. Comparison with sub-mm flux densities t owards mas- 
sive star formation regions (e.g. lHiU et aLll2005h . shows that 
in the worse case scenario, missing a 0.5 Jy contribution from 
free-free emission will result in an over-estimate of dust mass 
by a factor of ~2. 



5 CONCLUSIONS 

We have presented 19 to 93 GHz continuum and H70a -f 
H57q radio recombination line observations with the Aus- 
traha Telescope Compact Array towards 6 hot molecu- 
lar cores (HMCs) associated wit h methanol mas e r emis - 
sion. We confirm the results of iLongmore et al.l (|2007al ) 



that previous continuum surveys have missed several, some- 
times bright ('^O.SJy) c m-continuum sources. If the original 
iLongmore et al.l l|2007ah sample of 21 massive star formation 
regions are representative of the population as a whole, the 
evolutionary stage of a large number of regions have been 
mis-classifled. Rather than being very young objects prior to 
UCHii region formation, they are, in fact, associated with 
free- free emission and thus significantly older. 

In addition, inadvertently ignoring a potentially sub- 
stantial flux contribution (up to ~0.5 Jy) from free- free emis- 
sion has implications for dust masses derived from sub-mm 
flux densities, particularly the larger spatial scales probed by 
single-dish telescopes. Dust masses may be overestimated by 
a factor of --^2. 

We note that cm-continuum data from observations si- 
multaneously studying masers are particularly susceptible 
to this missing flux problem, due to the extended array con- 
figurations needed to achieve sufflcient angular resolution. 
The data in this work adds further evidence that high res- 
olution interferometric continuum data from such data may 
not be telling the whole picture. 



6 ACKNOWLEDGEMENTS 

We thank Paul Ho and the referee, Stuart Lumsden, for 
instructive comments. The Australia Telescope is funded by 
the Commonwealth of Australia for operation as a National 
Facility managed by CSIRO. 



REFERENCES 

Beltran M. T., Cesaroni R., Moscadelli L., Codella C, 

2007, A&A, 471, L13 
Cesaroni R., Hofner P., Walmsley C. M., Churchwell E., 

1998, A&A, 331, 709 
Forster J. R., CasweU J. L., 2000, ApJ, 530, 371 
HiU T., Burton M. G., Minier V., Thompson M. A., Walsh 

A. J., Hunt-Cunningham M., Garay G., 2005, MNRAS, 

363, 405 

Kurtz S., 2005, in Cesaroni R., FeUi M., ChurchweU E., 
Walmsley M., eds, lAU Symposium Hypercompact Hll 
regions, pp 111-119 

Kurtz S. E., Watson A. M., Hofner P., Otte B., 1999, ApJ, 
514, 232 

Longmore S. N., Burton M. G., Barnes P. J., Wong T., 
PurceU C. R., Ott J., 2007a, MNRAS, 379, 535 

Longmore S. N., Burton M. G., Barnes P. J., Wong T., 
PurceU C. R., Ott J., 2007b, in lAU Symposium Vol. 242 
of lAU Symposium, The molecular environment of mas- 
sive star forming cores associated with Class II methanol 
maser emission, pp 125-129 

Longmore S. N., Burton M. G., Minier V., Walsh A. J., 
2006, MNRAS, 369, 1196 

Longmore S. N., Burton M. G., PurceU C. R., Barnes 
P., Ott J., 2008, in Beuther H., Linz H., Henning T., 
eds. Astronomical Society of the Pacific Conference Series 
Vol. 387 of Astronomical Society of the Pacific Conference 
Series, Determining the Relative Evolutionary Stages of 
Very Young Massive Star Formation Regions, pp 58 — h 



Extended free-free emission towards MSF regions 9 



Longmore S. N., Maercker M., Ramstcdt S., Burton M. G., 

2007, MNRAS, 380, 1497 
Phillips C. J., Norris R. P., Ellingsen S. P., McCulloch 

P. M., 1998, MNRAS, 300, 1131 
Purcell C. R., Balasubramanyam R., Burton M. G., Walsh 

A. J., Minier V., Hunt-Cunningham M. R., Kedziora- 

Chudczer L. L., Longmore S. N., 27** other authors 2006, 

MNRAS, pp 198-+ 
Purcell C. R., Longmore S. N., Burton M. G., Walsh A. J., 

Minier V., Cunningham M. R., Balasubramanyam R., 

2009, MNRAS, 394, 323 
Sault R. J., Teuben P. J., Wright M. C. H., 1995, in Shaw 

R. A., Payne H. E., Hayes J. J. E., cds, ASP Conf. Ser. 

77: Astronomical Data Analysis Software and Systems IV 

A Retrospective View of MIRIAD. pp 433-+ 
Walsh A. J., Burton M. G., Hyland A. R., Robinson G., 

1998, MNRAS, 301, 640 
Walsh A. J., Macdonald G. H., Alvoy N. D. S., Burton 

M. G., Lee J.-K., 2003, A&A, 410, 597 



10 S.N.Longmore et al. 




Figure 1. G316.81-0.06 maps. The top-left image shows the integrated NH3(1,1) emission as dashed contours and the 24GHz continuum 
emission as sohd contours (taken from L07). The numbers show the labels given to each of the cores. The remaining images show 
the continuum emission at each frequency and array configuration. The peak flux shown at the top of each image corresponds to the 
maximum value in that image rather than the fits to the emission reported in Table [3] Contours are in 10% steps of the peak value with 
the maximum contour at 90%. The primary and synthesised beam sizes are shown as a circle and filled ellipse, respectively. The linear 
scale is illustrated by the line in each image. The position of the methanol maser emission and 8 GHz continuum from WBHR98 are 
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Figure 2. Radio recombination line spectra towards G316.81-0.06 in units of Jy/beam. Tlie spectra are extracted from the peak of the 
emission for each of the cores and are labelled by core number & transition in the top right, and array configuration in the top left of 
each spectrum. Overlayed are the Gaussian fits to each spectrum outlined in §|2]and used to obtain the values in Table |4] 
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Figure 3. G331.28-0.19 maps. The top-left image shows the integrated NH3(1,1) emission as dashed contours and the 24GHz continuum 
emission as sohd contours (taken from L07). The remaining images show the continuum emission at each frequency and array configura- 
tion. The peak flux shown at the top of each image corresponds to the maximum value in that image rather than the fits to the emission 
reported in Table|3] Contours are in 10% steps of the peak value with the maximum contour at 90%. The primary and synthesised beam 
sizes are shown as a circle and filled ellipse, respectively. The linear scale is illustrated by the line in the lower left of each image. The 
position of the methanol maser emission is shown as crosses, with the same angular size in each image. 
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Figure 4. Radio recombination line spectra towards G331.28-0.19 in units of Jy/beam. The spectra are extracted from the peak of 
the emission and are labelled by core number & transition in the top right, and array configuration in the top left of each spectrum. 
Overlayed are the Gaussian fits to each spectrum outlined in § [2] and used to obtain the values in Table |4] 
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Figure 5. G12.68-0.18 maps. The top-left image sliows ttie integrated NH3(1,1) emission as dashed contours and the 24GHz continuum 
emission as solid contours (taken from L07). The numbers show the labels given to each of the cores. The remaining images show 
the continuum emission at each frequency and array configuration. The peak flux shown at the top of each image corresponds to the 
maximum value in that image rather than the fits to the emission reported in Table [3] Contours are in 10% steps of the peak value with 
the maximum contour at 90%. The primary and synthesised beam sizes are shown as a circle and filled ellipse, respectively. The linear 
scale is illustrated by the line in the lower left of each image. The position of the methanol maser emission is shown as crosses, with the 
same angular size in each image. 
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Figure 6. Radio recombination line spectra towards G12.68-0.18 in units of Jy/beam. Tlie spectra are extracted from the peak of the 
emission and are labelled by core number & transition in the top right, and array configuration in the top left of each spectrum. Overlayed 
are the Gaussian fits to each spectrum outlined in § |2] and used to obtain the values in Table |4] 
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Figure 7. G19.47+0.17 maps. The top-left image shows the integrated NH3(1,1) emission as dashed contours and the 24GHz continuum 
emission as soHd contours (taken from L07). The numbers show the labels given to each of the cores. The remaining images show 
the continuum emission at each frequency and array configuration. The peak flux shown at the top of each image corresponds to the 
maximum value in that image rather than the fits to the emission reported in Table [3] Contours are in 10% steps of the peak value with 
the maximum contour at 90%. The primary and synthesised beam sizes are shown as a circle and filled ellipse, respectively. The linear 
scale is illustrated by the line in the lower left of each image. The position of the methanol maser emission is shown as a cross, with the 
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Figure 8. Radio recombination line spectra towards G19.47+0.17 in units of Jy/beam. The spectra are extracted from the peak of 
the emission and are labelled by core number & transition in the top right, and array configuration in the top left of each spectrum. 
Overlayed are the Gaussian fits to each spectrum outlined in § [2] and used to obtain the values in Table |4] 
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Figure 9. G24.79+0.08 maps. The top-left image shows the integrated NH3(1,1) emission as dashed contours and the 24GHz continuum 
emission as soHd contours (taken from L07). The numbers show the labels given to each of the cores. The remaining images show 
the continuum emission at each frequency and array configuration. The peak flux shown at the top of each image corresponds to the 
maximum value in that image rather than the fits to the emission reported in Table [3] Contours are in 10% steps of the peak value with 
the maximum contour at 90%. The primary and synthesised beam sizes are shown as a circle and filled ellipse, respectively. The linear 
scale is illustrated by the line in the lower left of each image. The positions of the methanol maser spots arc shown as crosses, with the 
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Figure 10. Radio recombination line spectra towards G24.79+0.08 in units of Jy/beam. The spectra are extracted from the peak of the 
emission for each of the cores and are labelled by core number & transition in the top right, and array configuration in the top left of 
each spectrum. Overlayed are the Gaussian fits to each spectrum outlined in §[2] and used to obtain the values in Table |4] 
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Figure 11. G24.85+0.09 maps. The top-left image sliows tlie integrated NH3(1,1) emission as dashed contours and the 24GHz continuum 
emission as solid contours (taken from L07). The numbers show the labels given to each of the cores. The remaining images show the 
continuum emission at each frequency and array configuration. The peak flux shown at the top of each image corresponds to the maximum 
value in that image rather than the fits to the emission reported in Table |3] Contours are in 10% steps (except the bottom three contours 
in the bottom right image which are separated by 5%) of the peak value with the maximum contour at 90%. The primary and synthesised 
beam sizes are shown as a circle and filled ellipse, respectively. The linear scale is illustrated by the line in the lower left of each image. 
The position of the methanol maser emission is shown as crosses, with the same angular size in each image. 
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Figure 12. Radio recombination line spectrum towards G24.85+0.09 in units of Jy/beam. The spectrum is extracted from the peak of 
the emission and labelled by core number & transition in the top right, and array configuration in the top left of the spectrum. Overlayed 
are the Gaussian fits to each spectrum outlined in § [2] and used to obtain the values in Table |4] 
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Figure 13. Spectral energy distributions for eacli of tlie sources detected at 24 GHz and not 8 GHz in WBHR98. Crosses and error bars 
show tile flux densities and tlreir corresponding uncertainties measured in tiiis work and L07A. To mitigate problems comparing fluxes 
of extended sources sampled at diflferent spatial frequencies, extended sources were re-imaged only including visibilities in the mutually 
overlapping spatial frequency range. Fluxes extracted from these images therefore sample the sources over similar spatial scales. The solid 
line shows the weighted best-fit to the detections and the corresponding spectral index is given as text. The arrows show the detection 
limit at 8 GHz reported in WBHR98. The dashed line shows the SED between 8 to 19 GHz of a source with a flux density at 8 GHz equal 
to the lower limit reported by WBHR98 and spectral index of 2. The horizontal dotted line shows the expected 19GHz flux that would be 
expected for a source of optically thick free-free emission, which lies just below the detection threshold of WBHR98. The 8 GHz detection 
limit and dashed line are omitted for G19. 47-1-0. 17 as this was too far outside the WBHR98 primary beam to determine a reliable 
detection limit. The dashed lines in the remaining sources show that, with the exception of G24.79+0.08, the measured flux densities in 
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